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Abstract 

We show that the leading power corrections to the event shape distributions can be 
resummed into nonperturbative shape functions that do not depend on the center-of-mass 
energy and measure the energy flow in the final state. In the case of the thrust variable, 
the distribution is given by the convolution of the perturbative spectrum with the shape 
function. Choosing the simplest ansatz for the shape function we find that our predictions 
for the thrust distribution provide a good description of the data in a wide range of energies. 
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1. Introduction 



Study of hadronization corrections to the event shapes in e"*" e~ annihilation became a unique 
laboratory for testing QCD dynamics beyond perturbative level . Being infrared and collinear 
safe quantities, the event shapes (their mean values as well as differential distributions) can 
be calculated in perturbative QCD at large center-of-mass energies s = Q"^ as series in as{Q). 
Nonperturbative corrections to the event shapes are attributed to hadronization effects and they 
are expecting to modify perturbative predictions by terms suppressed by powers of large scale 
with the exponent p varying for different observables. 

Successful description of the hadronization effects by the phenomenological Monte-Carlo based 
models indicates that in distinction with the total cross-section of e"*" e~ annihilation the power 
corrections to the event shapes become anomalously large and for the shape variables like the 
thrust and heavy mass jet they are expecting to appear at the level p = 1. 

The enhancement of hadronization corrections occurs due to the fact that the event shapes 
are not completely inclusive quantities with respect to the final states but rather weighted cross- 
sections in which large power corrections can be attributed to an incomplete cancellation of the 
contribution of soft gluons. As a consequence, the operator product expansion (OPE) is not 
applicable to the analysis of the event shapes and the standard identification of the exponents 
p characterizing the strength of power corrections as dimensions of local composite operators 
entering the OPE does not hold. 

To determine the leading exponent p and also to understand the way in which nonperturba- 
tive effects modify perturbative predictions one may explore instead by now standard infrared 
renormalon analysis. This procedure has been successfully applied to the mean value of different 
event shapes variables and the description of the leading power corrections has been given 
within different approaches 0-0. In contrast, the hadronization corrections to the differential 
event shape distributions are less understood. One of the reasons for this is that the leading 
power corrections to the mean values and to the differential distributions have different form 
the former are characterized by a single nonperturbative scale of dimension p while the latter 
involve the nonperturbative function of the shape variable that one usually estimates running 
the Monte-Carlo event generators. 

Studying the power corrections to the event shape distributions we will follow the approach 
proposed in 0. We will mostly concentrate on the differential distribution with respect to the 
thrust variable da/dT and, particularly, in the end-point part of the spectrum T ~ 1. 

There are few reasons for considering the region T ~ 1. In contrast with the mean value 
(1 — T) that gets l/Q— power correction from the final states with an arbitrary number of jets, 
for the thrust distribution in the end-point region, T — 1, one has in the final state only two 
narrow energetic jets moving close to the light-cone directions p+ and p_ {p± = ^(1, 0, ±1) and 
2(p+p_) = Q"^) in two opposite hemispheres. Denoting their invariant masses as M| and one 
gets 



where for later convenience we introduced a new variable t. 

Taking into account the QCD effects of collinear splitting of quark and gluons inside two 
narrow jets and their interaction with soft gluon radiation one finds that the thrust distribution 
for t — i> depends on two infrared scales, Q'^t^ and Q^t, which give rise to large both perturbative 
(Sudakov) logs and power corrections. The smallest scale, {QtY, sets up the total energy carried 
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by soft gluons in the final state and the scale Q'^t characterizes the transverse size of the jets 
/c^ = 0{Q'^t). The power corrections to the thrust distribution are suppressed by powers of both 
scales. In order to separate the leading asymptotics one keeps the smaller scale Qt fixed and 
expands the thrust distribution in powers of the larger scale Q^t 

^ '^'^ aJa,{Q)Mt,^\+o{^\. (2) 



cjtot dt \ ' ' Qt) \Q^t^ 

In what follows we will consider only the leading term of this expansion, ctq. It should resums 
large perturbative terms Oi^{Q) In"^ t/t {m < 2n — 1) and take into account all power corrections 
of the form l/iQi)"^. 

The structure of power corrections ~ l/lQt)"^ strongly depends on the value of the thrust 
variable. Away from the end-point region, t ^ A^^^/Q one may retain in ctq only the lowest term 
n = 1 and neglect the terms with n > 2 as suppressed by powers of l/Q- It is this approximation 
that one applies calculating the mean value of the thrust (t) = a^^t /o™'"' dtt^ , where tmax is the 
upper limit on the thrust variable that one imposes to separate the contribution of the 2-jet final 
state, tmax = |- At the same time, in the end-point region t = 0{A^^^/Q) all terms ~ {Qt)~^ 
become equally important and need to be resummed inside ctq for all n. 

As we will show, the leading power corrections to the thrust distribution in the end-point 
region away from the small invariant jet mass limit Q^t ^ ^qcd '^^^ resummed into a non- 
perturbative Q— independent function that defines the shape of the distribution in the region 
t = 0{A^^^/Q) and therefore is called the shape function. Then, the QCD prediction for the 
leading term in (0) is given by the convolution of perturbative Sudakov spectrum with nonper- 
turbative shape function. One should mention that one finds similar expressions considering, for 
example, the large— x asymptotics of the structure function of deep inelastic scattering [Rl M and 



the end-point spectrum of the inclusive heavy meson decays B 7X5. |T0|, [TT|, [T^. The reason 
for this similarity is that in all these cases one encounters the same physical situation when en- 
ergetic narrow jet(s) is propagating in the final state through the cloud of soft gluons. However 
the important difference with the thrust distribution is that the nonperturbative functions in 
the latter two cases resum power corrections on the different scale Q^t and they coincide with 
well-known inclusive (light-cone) distributions. 



2. Analysis of soft gluon effects 



At the Born level the final state consists of a quark-antiquark pair and the thrust distribution 
has the form 

—-J7 =m- (3) 

O"tot Born 

Soft gluon radiation smeares this peak towards larger t. Let us first analyze separately pertur- 
bative and nonperturbative contributions. 

Considering the perturbative emissions of soft gluons out of two outgoing quarks one finds 
that for t — i> the phase space for real soft gluons is squeezed and due to an incomplete can- 
cellation between virtual and real gluon contributions the perturbative corrections to the thrust 
distribution involve large Sudakov logs that can be resummed to all orders with the double 
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logarithmic (DL) accuracy as 
1 da 



(Ttot dt 



PT 



i?,,(t)=exp In t (4) 



with Rprj,{t) called the radiator function. One can systematically improve perturbative approxi- 
mation by including additional nonleading logarithmic terms in /2p^(t) and matching the result 



into exact higher order calculations using the In i?— scheme [|T3|. The perturbative Sudakov spec- 
trum extends over the interval < t < tmax and vanishes at the end points. The peak of the 
distribution is located close to t = and it is shifted towards larger t as one improves pertur- 
bative approximation. Its position, tp = 0{A^^^^/Q), is sensitive to the emission of soft gluons 
with energy ~ A^^.^ indicating that the physical spectrum around the peak is of nonperturbative 
origin. 

Let us now estimate the effects of nonperturbative soft gluon emissions on the thrust dis- 
tribution @. We take into account that in the leading order in l/{Q'^t) the transverse size of 
two quark jets k'^ = 0{QH) can be neglected, that is soft gluons with the energy ~ Qt can not 
resolve the internal structure of jets. This means that considering soft gluon emissions we may 
apply the eikonal approximation and effectively replace quark jets by two relativistic classical 
particles that carry the color charges of quarks and move apart along the light-cone directions p+ 
and The interaction of the quark jets with soft gluons is factorized into the unitary eikonal 
phase Wifi) given by the product of two Wilson lines calculated along classical trajectories of 
two particles 

Vr(0) = ly+(0)[l^_(0)]^ W±{x)= Pe^^[i dsp±- A{x + p±s)^ , (5) 

with gauge fields A^{x) describing soft gluons. Denoting the total momentum of soft gluons 
emitted into the right and left hemispheres as kji = J2ieR^i and ki = J2ieL^i: correspondingly, 
one finds the thrust (0) as t = 2{kRp+)/Q'^ + 2{kLP-)/Q^ and obtains the following expression 
for the differential distribution 

' j:\{mmN)f5(t-^-%] (6) 



c^tot dt N \ Q Q 

with fc^ = ko^k^. Here, the matrix element of the Wilson line operator describes the interaction 
of quarks with soft gluons and the summation goes over the final states N of soft gluons with 
the total momentum k = k^ + ki. Expression (P) follows from the universality of soft gluon 
radiation and it takes into account both perturbative and nonperturbative corrections p. 

Let us neglect for the moment the perturbative contribution to the matrix element of the 
Wilson line in (^. Then, introducing the shape function 

m = j:mwmN)\"s{e-k],-kt) (?) 

N 

one can estimate the nonperturbative contribution to the thrust distribution as 

— ^ =Qfm- (8) 

O^tot nonPT 

The nonperturbative function f{e) is localized at small energies e and according to ® it deter- 
mines the shape of the spectrum at small t = 0{A^^^/Q). 
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The important property of the function (|^) is that it does not depend on the center-of-mass 
energy Q. Although the vectors p± entering into the definition (|^) of W(Qi) do depend on Q, 
the Q— dependence of the Wilson lines (^) disappears due to the reparameterization invariance 
s —>■ As. Therefore, one may extract the shape function at some reference Qq and then apply it 
to describe the hadronization effects for different center-of-mass energy. 

Being a new nonperturbative distribution, the shape function admits the operator defini- 
tion that is different however from the one for the inclusive distributions. As a manifestation of 
noninclusiveness of the thrust variable, f{e) depends separately on soft gluon momenta flowing 
into different hemispheres. In particular, it takes into account the effects of soft gluon splittings 
when the decay products fly into two different hemispheres 0. Since this may happen at time 
scales larger then 1/Q, we do not expect the function f{e) (more precisely its moments) to be 
related to matrix element of local operators as it happens for inclusive distributions. Indeed, the 
operator definition of /(e) involves the "maximally nonlocal operator" S{n) that measures the 
density of the energy flow in the direction of unit 3- vector n. According to its definition S{n) acts 
on the final state of particles as £{n)\N) = Y^iLi k^5{n — ki/k^)\N) and it can be expressed in 
terms of the tensor energy-momentum operator [0, Then, one uses to get 

f{e) = {0\W^{0)6(^e- Jd'n{l-\cos9ri\)S{n)^W{0m, (9) 

where 9fi is the angle between vector n and the thrust axis, cos 9fi = n^. The detailed properties 
of this function will be discussed elsewhere WB . 



3. Factorization of the thrust distribution 

The expression (||) was found by neglecting perturbative soft gluons. It is clear that they also 
contribute to the matrix element of Wilson lines entering and modify the thrust distribution. 
In order to combine together perturbative and nonperturbative effects one has to introduce the 
factorization scale n and separate the contribution of soft gluons with the momentum above and 
below this scale into perturbative Sudakov spectrum, and nonperturbative shape function, 

/(e), respectively. Both quantities become functions of the IR cut-off n but the thrust distribution 
is /i— independent. 

For the inclusive distributions the above procedure can be performed using the operator 
product expansion. For the thrust distribution the OPE is not valid and we apply instead the 
infrared renormalon analysis. To this end we perform perturbative calculation of (|^) by summing 
over the final states | A^) of multiple perturbative soft gluon radiation and identify the ambiguities 
of resummed perturbative expressions that can be attributed to nonperturbative contribution. 
One finds that thanks to the nonabelian exponentiation of the Wilson lines, the perturbative soft 
gluon contribution exponentiates in the Laplace transform of the distribution |Q 

(e--*) ^ / dte-'^' —— = exp {-S{u)) (10) 

Jo (Ttot dt 

with the leading term in the exponent of the following form 



5 



and rcusp(as) = ^ + C'lttg) being a universal cusp anomalous dimension. 

We observe that since the singularities of the coupling constant affect the integration over 
transverse momenta of soft gluons in ( pJ]) the Sudakov form factor ^(z/) suffers from infrared 
renormalon ambiguities. They originate from soft gluons with the energy of order Aqcd whose 
contribution should be separated into nonperturbative function f{e). Namely, introducing the 
cut-off /i on the value of transverse gluon momenta in ([TTD we may split S{i>) into the sum of two 
terms. The term with > /i^ defines the perturbative contribution to the exponent, S'p^(i/), 
which in turn allows to find the perturbative spectrum '^'^pt^^'^^^ through the inverse Laplace 
transformation (|10|). The second term with k'j_ < /i^ should be absorbed into the definition of the 
nonperturbative function (0). In this case, changing the order of integration in ([TT|) one expands 
the integral in powers of z/ and absorbs the ambiguous integrals Jq dk±k']^^Tcnsp{c(s{k±)) into the 
definition of new nonperturbative scales An(/i). Substituting (H) into ( p!0D one finds the following 
consistency condition 

/ c?ee— = exp - J] (z//g)"A„(/x) + 0{u/Q') . (12) 

•^0 \ n=l ) 

Although we can not draw any conclusions about the absolute value of the scales A„, their 
/i— dependence is of perturbative origin and it can be determined as 

^dX^ = 4(-)"+in- V'^r.^splasl/i)) . (13) 



Since the parameter v is conjugated to the thrust t we neglected in (|T^ the corrections ~ vjQ^ 
and replaced the upper limit of e— integration, ^max = ^ma^^Q-, by oo. 

The fact that the infrared renormalons contribute additively to the exponent of (^) implies 
that the Laplace transform of the thrust distribution is factorized into the product of perturbative 
and nonperturbative terms 

/■oo 

(e-*) = (e-*)^^ X / d^e— /«/(e;/.) (14) 

JO 

where (e~'^*)p^ is calculated as a mean value with respect to the perturbative distribution 
Integrating the both sides of this relation with respect to v with an appropriate weight we obtain 
the factorized expression for the radiator function -R(t) = /q^ dt^ = {9{t — t)) 

R{t) = de f{e; fi)R^^ (^t " |; , (15) 

where the upper limit of integration follows from the condition i?p^(t) = for t < 0. Thus, the 
net effect of incorporating nonperturbative corrections (in the leading 1/ (Q'^t)— order) amounts 
to the l/Q— shift of perturbative radiator function smeared with the shape function. 

To see how (|15|) resums both perturbative and nonperturbate corrections, one expands the 
radiator R-prj,(t — ^] fdj in powers oil/Q 

R{t) = R,,it) - Mi?;^(t) + ^[K't(^) - KM + • • • ' (16) 
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where prime denotes the logarithmic derivative with respect to t. Here, the leading term Rprj,{t) 
gives the resummed perturbative Sudakov expression which is different from due to additional 
dependence on the IR cut-off fi. The terms with the derivatives of Rp^ generate the series of power 
corrections accompanied by the set of nonperturbative /i— dependent dimensionful parameters 
(e") that can be expressed in terms of the scales introduced in (|T2|) as 



(e")=/ {e) = , {e)' - {e') = X2 , ... (17) 



Finally, differentiating the both sides of (|r^) with respect to t and taking into account that 
Rp^{t] /i) = 9{t) Jq (lt^lfp:^2tl find the thrust distribution 

1 da(t) /"Q* dapJt - f^; a) 

^ = Q/(Qt;/iR,(0;/i)+ / de^\e■^^)^^—-^^. (18) 

a tot dt Jo dt 

Here, the two terms in the r.h.s. have the following interpretation. The first term dominates 
at small t and corresponds to the situation when the perturbative real soft gluon radiation is 
washed out due to the IR cut-off and the final state consists entirely of nonperturbative radiation 
described by the shape function. In comparison with (H) one gets the additional Sudakov factor 
Rprp{0; fi) that takes into account the perturbative contribution of virtual soft gluons with the 
energy above the cut-off fi. This factor rapidly vanishes as one decreases the value of fi allowing 
more perturbative virtual gluons to be emitted. The second term in (|18D describes the smearing of 
the perturbative Sudakov spectrum by nonperturbative corrections. As an illustration of (0), we 
depicted on Fig. 1 the contribution of both terms to the thrust distribution at the center-of-mass 
energy Q = 35 GeV. 

In contrast with the heavy meson decay where nonperturbative corrections extend the per- 
turbative spectrum beyond the perturbative end-point due to interaction of the heavy quark in 
the initial state with the light component of the meson |jTO|, |ll|, the nonperturbative corrections 
to the thrust distribution have just an opposite effect. They shift the perturbative spectrum 
inside the perturbative window < t < tmax and describe an "evaporation" of the energetic jets 
in the final state due to emission of soft gluons. 



4. Nonperturbative ansatz for the shape function 

The shape function is a new distribution that resums nonperturbative corrections to the thrust 
distribution. Although its explicit form can not be extracted from our analysis we could use the 
renormalon inspired sum rules (|1^) to study its general properties. 

Expanding the both sides of ([I^ ) in powers of 1/Q we verify that in accordance with its 
operator definition (^ the function f{e) does not depend on the large scale Q and its first few 
moments are given by (plT]) . The shape function depends on the factorization scale fi and one 
may apply the renormalization group equations (|T3|) to find its evolution with /x. 

Using ([T7|) one may formally write the shape function in the form of the distribution as a 
series in 5— function and its derivatives 

f{e; ft) = 5{e - Ai) + ■ 5'ie - X,) - h^d'^e - X,) + ... (19) 

Its substitution into (|T8|) yields a series that is equivalent to the expansion of the radiator (|I6|) . 
It is convergent however only for the values of the thrust A„^^/Q <^ t < tmax on the tail of the 
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NLL+power corr. 

NLL 




Figure 1: The prediction (|T8D for the thrust distribution at Q = 35 GeV. The dotted and dashed 
hues describe the first and the second term in the r.h.s. of (|T^), respectively, and the sohd hne is 
the sum of both. The dot-dashed hne denotes the perturbative Sudakov spectrum 
the NLL approximation. 



do-prp (t;/j=0) 



m 



Sudakov spectrum where the perturbative distribution is a slowly varying function of t. In 
this range of t, keeping only the first term in the r.h.s. of ([T9|) , one finds that the leading power 
corrections simply renormalize the thrust variable by generating the shift of the perturbative 
spectrum [Q 



1 da{t) da^^{t-^' 



(Ttot at dt 

This result is a general property of the leading l/Q— power corrections to different event shapes 
and it is an immediate consequence of the exponentiation of soft gluon contribution. The pre- 
diction ( pOD has been found to be in a good agreement with the data [|1^. 

For the values of the thrust t = 0{Af^f^j^/Q) one needs to know the explicit form of the shape 
function. In this case, one relies on the particular ansatz for f{e) that can be inspired by different 
model considerations. As the simplest form of f{e) we choose the following one 



fie) 



exp 



. AVAr( 



(21) 



A 



This function depends on two parameters: dimensionless exponent a controlling how fast the 
function vanishes at the origin and dimensionfuU scale A defining the interval of energies on 
which the function is localized. The shape function is peaked around ^max 
rapidly vanishes for e > ^max- 

There are additional constraints that we may impose on the shape function. They follow 
from the analysis of nonperturbative corrections to the mean value of the thrust and its higher 
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moments. Indeed, expanding the both sides of ( p!4D in powers of v one gets 

(t) = (t)p. + ^, (t^) - (t)^ = - (t)^. + ^^^^^ , - (22) 

These relations are vahd up to corrections due to contribution to the thrust of the final states 
with n > 3 jets. Using (p2D and ([T7| ) we may relate the parameters of the shape function to the 
mean value of the first few moments of the thrust. In particular, applying the results of ^ on 
parameterization of l/Q— corrections to the mean value of the thrust we get 



1 



0.8 GeV 



Qh ' Q 

Substituting ( pT]) into this relation one still has a freedom in choosing the single parameter a 
that we define as 

a = 3, A = 0.7 GeV. (23) 

The shape function corresponding to these values of parameters is shown in Fig. 2. One should 
stress that the explicit form the shape function ( pi] ) as well as the values of the parameters a 
and A are related to the particular choice of the factorization scale /x to be specified later on and 
they do not depend on the center-of-mass energy Q. 




e (GeV) 



Figure 2: The nonperturbative ansatz for shape function. 

Let us now consider the perturbative spectrum ^^^pii^i/fl^ 'jj^g contribution of perturbative soft 
gluons with the energy above the IR cut-off /i exponentiates in the Laplace transform ([T0|) and 
S'p^(z/) is given by (|ll]) with the additional condition k\ > /i^ imposed on the fc_L— integration. 
The important difference with the known results for resummed Sudakov spectrum [0] is in the 



additional /z— dependence of ^fp:^2t}._ Nevertheless, one may expand following the exponent 
5'p^(z/) in powers of a^iQ) and separate large logarithmic terms a^U^ {m < n + 1) which due 
the presence of additional scale appear of two kinds, L = hi[Q/ii) and L = ln(l/z/). Then, in 
the NLL approximation that takes into account all terms a^L"""*"^ and a"L" in S'p^(i/) one may 
replace 1 — e~"'^ ^= 6{u — e~'^E /i/) in (|l^) and find that due to the condition k±> fi the function 
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Sprj^lu) has different behaviour depending on the value of u, or equivalently the thrust t. Finally, 
one performs the inverse Laplace transformation and obtains the radiator function in the NLL 
approximation as 

p ,A NLL exp(-gp^(to/t)) dSp^jto/t) 
^) = T{l-S'^^{to/t)) ' = d\nt ' 

with to = ^"'^ = 1.780.... Examining this expression one finds that for t/tQ > fi/Q the radiator 
Rp^{t; n) rapidly vanishes at small t and it does not depend on the IR cut-off /i. Thefefore for 



these values of t it coincides with the known expression from m. For /^VQ^ <.t/tQ< ^i/Q the 
radiator Rp^it) starts to depend on the IR cut-off yU and, as a consequence, its decrease at small 
t is slowed down. For < t/tQ < fi'^/Q'^ the radiator is t— independent. 

To simplify numerical calculations of the spectrum (|1^) we ignore the difference between 
Rprj,(t = tofi/Q) and Rprj,(t = t^fi'^/Q'^) by choosing the factorization scale /i to be small but 
within the applicability range of the NLL approximation [|1^], 2Poas{Q)ln{Q/ fj,to) < 1. In this 
case, one may approximate the radiator (|2^) as 



RZ^it; fi) = 9{t/to - fi/Q)Rp^{t) + 9{t)9{fi/Q - t/to)Rp^{fi/Q) , (25) 



where i?p^(t) is the known expression for the radiator in the NLL approximation improved 
by higher order corrections in the In i?— matching scheme. We choose the two free parameters, 
the factorization scale and the fundamental QCD scale as 

^ = 0.750 GeV , Aqcd;ms = 0-250 GeV . 



Differentiating the expression ( psf ) with respect to t we obtain the perturbative spectrum 



that starts at t = t^^/Q and extends to t = t„ 



dt 



One should notice that if instead of (^) we would have used for t/t^ < /i/Q the expression 



for the radiator Rp^{t; fi) in the NLL order, (|2^, then the perturbative spectrum will have 
discontinuities at t/to = n/Q and t/t^ = /i^/Q^. The reason for this is that ^"'pT^^'f^^ = ^p^^^'^^') 
involves the second order derivative of the Sudakov form factor, S'^^, that, as it can be seen from 
(pTD, is not enhanced by large logarithmic corrections and therefore can not be approximated by 
the NLL result in the transition regions around these two points. 

Finally, let us compare the QCD prediction ([T8| ) with available data on the thrust distribution 
in the interval < t < 0.33 at different energies 14 < Q/GeV < 161. Fig. 3 shows the comparison 
with the data at Q = 91.2 GeV, where the most accurate experimental data are available. The 
combined fit for various center-of-mass energies |T3] is shown in Fig. 4. We would like to stress 



that for different values of Q we use the same ansatz for the shape function, (|21|), with the 
parameters defined in (|23|) . As a nontrivial test of ([T8|) we observe that the theoretical curves 
reproduce well the Q— dependence of the data in the end-point region. 



5. Conclusions 

We have study the power corrections to the thrust differential distribution in the end-point region 
of t. Due to enhancement of soft gluon contribution, the spectrum is affected by large perturbative 
Sudakov and nonperturbative l/(Qt)— corrections that need to be resummed. Using universal 
properties of soft gluon radiation we have argued that the resummed leading power corrections 
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Figure 3: The comparison of the data with the QCD prediction for the thrust distribution ([18|) 
at Q = 91.2 GeV 

are described by the shape function which is a new nonperturbative distribution independent on 
the center-of-mass energy Q measuring the energy flow in the final state. The thrust distribution 
is given by the convolution of the shape function and the perturbative Sudakov spectrum each 
depending on the factorization scale /x. For the values of the thrust t > ^/Q the power corrections 
generate the shift of the perturbative spectrum, (pOl) . In this case, the thrust distribution is not 
sensitive to the particular form of the shape function but only to its first moment (e). In contrast, 
for A^^^/Q^ ^ t < yu/Q it is the shape function that governs the end-point behaviour of the 
spectrum. Choosing the simplest ansatz (^) for this function and using the In i?— matched 
expression for the perturbative spectrum ( ^5]) we have found that our prediction (0) provides a 
good description of the data in a wide range of energies. 

Analysing the power corrections to the event shapes one should identify universal nonpertur- 
bative quantities that describe the hadronization effects in e"*" e~ annihilation. According to (|^) 
the shape function depends on the definition of the thrust variable and considering the power 
corrections to other event shapes like heavy mass jet or energy-energy correlations one gets the 
expressions for the shape functions as well as the factorized expressions for the distributions dif- 
ferent from (P) and ([T8|) . However, taking the moments of the shape function, / dee^ f{e), one 
finds that for various event shapes they are expressed in terms of the same universal distribution 
{Q\W\^)E{ni)...£{nN)W{<d)\^) that measures the energy flow in the final state in the directions 
specified by unit 3- vectors wi, n^. This object deserves additional studies [|T6]. 
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Figure 4: The comparison of the data with the QCD prediction for the thrust distribution at 
different energies (from bottom to top): Q/GeV = 14 , 22 , 35 , 44 , 55 , 91 , 133 , 161. P| 
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